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Abstract 

It is shown that the Kapitza-Dirac effect with atoms, which has been 
considered to be evidence for their wavelike character, can be interpreted as 
a scattering of pointlike objects by the periodic laser field. 

1 Introduction 

The currently accepted answer to the question posed in my title is, of course, 
"Both". But I submit that we should not abandon the heritage passed down 
to us by the Atomists, from Democritus to Boltzmann. It was a long struggle, 
at times involving scientific isolation and consequent personal suffering^], to 
establish the atomicity of matter (from which I exclude radiation for present 
purposes) . 

In recent times some experimental evidence has been found to support 
a wavelike description of atoms, and even of quite large molecules such as 
fuUerene. I shall concentrate here on the first category, in which something 
analogous to the diffraction of light has been observed with a "monochro- 
matic" beam of atoms, the grating being supplied by a stationary laser wave 
which is tuned to a frequency close to an atomic resonancej21. Actually it 
is not so much the monochrome property which is important - the velocity 
of the beam was controlled only to within about 5% of its mean value - but 
rather a very high degree of collimation - the component of momentum per- 
pendicular to the laser must be an order of magnitude less than h/X, where h 
is Planck's constant and A is the laser wave length. For sodium atoms with a 
mean velocity of lO^ms"^, and with the laser tuned to the D-line at 589nm, 
this demands an initial angular spread less than 3.10~^radians. What we 
observe in the outgoing beam is a set of well separated peaks at integral mul- 
tiples of 6.10~^rad. There is a well worked out theory of the broadening of 
these linesj2]|l], but the spacing of 2h/XMv, as well as the intensities, may 
be explained with a very simple quantum mechanical (QM) model to be out- 
lined in the next section. This model was discussed by Gould who offered 
two interpretations of the analysis; we must accept either that each atom is 
spread out over a wave front of the order of several microns or that the atom 
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trades in quanta of momentum. The first alternative is the description given 
long ago by Kapitza and DiracfS^ of a diffraction process in which an atomic 
wave whose wavelength is h/Mv is diffracted by a grating whose spacing is 
A/2, while the second views the process as one of scattering in which the 
atom absorbs and emits, stochastically, radiation from and to the laser field 
in quanta of he/ A; such radiation is in one of the two (up or down) directions 
of the laser beam, and so its Poynting vector carries a transverse momentum 
of ±h/X, and this must be compared with a longitudinal momentum of Mv. 
But, curiously, the latter analysis indicates that events of emission and ab- 
sorption occur in pairs, so that very few atoms emerge from the laser having 
a transverse momentum which is an odd multiple of h/X. 

Gould did not choose to emphasize the contradictory nature of these two 
interpretations; he instead pointed out that either of them were "equally un- 
palatable to the prequantum physicist" . Staying within the Atomist tradition 

1 propose to reject the first interpretation and accept the second. Neverthe- 
less, I enter the reservation that I can do so staying largely within a classical 
(or prequantum) world view. There is a fair amount of evidence jU] that Max 
Planck, who discovered the quantum discontinuity in absorption and emis- 
sion of light, never accepted that the light field itself had to be quantized. 
A quotation from a letter to Einstein in 1907 illustrates Planck's view of the 
light field. 

I am not seeking the meaning of the quantum of action (light 
quantum) in the vacuum but rather in places where emission and 
absorption occur, and I assume that what happens in the vacuum 
is rigorously described by Maxwell's equations. 

In summary, I propose, from Section 3 onwards, to investigate whether the 
distribution pattern of the scattered atoms may be explained on the basis of a 
model in which quanta h/X of momentum are exchanged, stochastically, with 
the laser field. Before that I shall summarize, in Section 2, the results of the 
simplified QM theory, which will provide us with a standard for comparison. 

2 The QM model 

The hamiltonian is 

H (t) = —h{uj — A) o"3 + hOsB, cos C ((Ti cos ujt + cr2 sin ut) , < t < to , (1) 
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where uo is the laser frequency, detuned by A from the D-hne resonance Uq, 
is the resonant Rabi frequency of the interaction, and cti,(T2,o'3 are the 
Pauh spin matrices interpreted, in the standard manner, as atomic raising 
and lowering operators. The kinetic energy [h? /2M}?) [d"^ / dC,"^) has been 
discarded on account of the large atomic mass M, and the variable C, which 
takes values in the range (— 7r/2,7r/2) (see Fig.l), gives the phase of the 



Figure 1: The Kapitza-Dirac effect according to the QM model, which is 
the original description of Dirac and Kapitza. A de Broglie wave, whose 
wavelength is small, and whose coherence width is large, compared with 
A, is incident on a stationary laser of wavelength A. The de Broghe wave 
experiences the laser as a sinusoidally varying refractive index with spacing 
A/2, and emerges as a diffraction pattern whose maxima are separated by 
angular intervals of 2h/XMv radians. 

atom in the laser wave at the point of entry. On account of the assumption 
of infinite mass, this is also the phase at the point of exit. Starting from an 
initial state having zero transverse momentum and in the lower state of the 
D-line couplet, that is 



V 



X 





(2) 



the state at time to is 



( 



[-i{nR/n) COS C sin fito] 6^'^*°^^ 
[cos nto + i (A/20) sin Qto] 6'^^°/^ 
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where = y^J^J^)~+T^co^^. The amphtudes '0u and ^pi are the parts of 
the wave function representing an atom in its upper and lower state respec- 
tively. In order to observe the discrete momentum spectrum it was found 
necessary to make A substantially larger than flu; typically A ?a 5^2^^. We 
therefore expand the wave function in powers of 7 = fin/ A giving, to order 

7^ 

i^uito'X) = — 2^7 cos C sin f2'r , 

(to; C) = e*^'" - 2^7^ cos^ C sin Q'r , (4) 

where 

T = WA/2 , = 7-2 + 1 + cos 2C . (5) 

The Fourier series for the transition amplitude is then 

(r; C) = -il J2 t"^- + ■^"+1 ^^'"^'^'^ ' (6) 

n=— c« 

where 

(l + 72)r^n7r 



and for the no-transition amplitude it is 



r 2 ' 



i^i (r; C) = E e^^"^ [e'^"Jn (r) - i7'{ J„ (r) sinr^ - 4 (r) cosrj] . (8) 
—00 

The lower component gives the intensities of the even lines of the spectrum, 
namely 

(r) = Jn ir) (1 - 272 sin^ r„) + 7V, (r) j; (r) sin 2r„ , (9) 
while the upper component gives the odd lines, namely 

Pn+1/2 i^) = 7^ [Jn (t) sin Tn + Jn+1 (t) cosr„]^ . (10) 

The wave interpretation of Kapitza and Dirac is made by considering the 
limit 7 — > 0, so that the outgoing wave function is effectively the scalar 



^ (^0; C) = exp 



, 1 2 > 

i [ —-ujoto + 2rcos C 
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The dependence oi ip on ( is an indication (see Fig.l) that the de Broghe 
wave of the atom experiences a spatially varying refractive index as it goes 
through the laser, and its Fourier expansion, that is 

oo 

exp {2tT cos' C)= Yl ^'"^ ' (12) 

ra=— oo 

indicates that the atom acquires a transverse momentum of either 2nh/X or 
—2nh/X with probability 

(r) = Jl (r) , (13) 
for which the characteristic function is 



Fqo {9- t) = (e^"^> =Y.Jl (r) e^' = J, [2r sin 



(14) 



The moments of the distribution are obtained from the derivatives of F. 
In particular the variance is 

oo ^ 

(^'>Q0 = E """Jl ir) = -FH,, (0; r) = -r' . (15) 

ra=— oo 

This has the form, for small r, 

(n'>Q„ = 2p?°(r)+0(r^) , (16) 

which establishes that changes in n occur in single steps of ±1. 

We shall need to consider the asymptotic behaviour as r ^ oo of this 
spectrum, namely |H] ^ 



27r-i (r^ - n')-^/' cos' [Vr^ - - P\n\ - 7r/4] (|n| < r) , 
(|n| > r) , 

(17) 

where 

/5 = cos-i (|n|/r) . (18) 

These intensities oscillate, with angular frequency 1/2 for small n, but de- 
creasing as n approaches r. However, the oscillations disappear once we take 

"'^In the transition region n ^ t these asymptotic expressions should be replaced by 
others, obtained from Airy approximations and also given in [H]. However, the velocity 
averaging which I propose next will mask this correction. 
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Figure 2: The momentum spectrum, averaged over the velocity profile, of 
the Kapitza-Dirac effect according to the QM model, with the time param- 
eter r = 50. The bar chart represents the intensity of the nth line and the 
continuous curve depicts a deterministic classical model. Since the observed 
datum is actually the angular deflection, the experimental method used can- 
not distinguish the two spectra at this value of r. 



account of the beam's velocity profile, described by a gaussian function H{t), 
with standard deviation a =0.025r, so that 95% of the atoms have transit 
times within 5% of the mean. The averaged intensities are 



r Jn ir') H (r, r') dr' \= [ 

Jo Txay iTX Jn 



exp 



dr' 



■/2 



_ (19) 
I have plotted, in Fig.2, the values of and its asymptotic limit at r = 50. 
The latter curve actually coincides with a completely classical model of the 
process, in which an atom going through the laser at phase C acquires a 
momentum of rsin2C, and C, has uniform density between —7^/2 and 7r/2. 
For such large r it is this classical curve which would be observed, because 
the experimental datum is the angular deflection, rather than the transverse 
momentum, of the atom; the individual lines of the spectrum are broadened, 
so that they merge with one another. 
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The exact spectrum, given by © and (fTUI). displays rapid oscillations 
because of the sinusoidal terms in r„, in addition to the slower oscillations 
we have just been considering. But, except in a short initial period, the 
rapid oscillations disappear for all n after smoothing with H. I have plotted 




0.5 1 1.5 



Figure 3: The intensities of the first few lines, averaged over the velocity 
profile, in the QM model for small values of the transit time r. The parameter 
7 has been set at 0.2. 

the smoothed spectrum for the first few values of n in Fig.3; with 7 = 
0.2 the rapid oscillations are effectively damped out for r > 1. For the 
actual experimental range of 2 < r < 6, we may smooth by putting simply 
(2sin^7:„) = (2cos^r„) = 1 and (sin2r„) = 0, leading to 




and this corresponds to the characteristic function 



1+7 



COS--I 



(21) 
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3 A stochastic model 



In his article, GoulcH2j states that "... if we attempt to assign specific points 
in the diffraction pattern to specific locations in the standing wave, we will 
fail miserably". While not dissenting from this judgement, I stress that 
many areas of classical physics produce situations of the same character; 
if we were to try to predict the position of a Brownian particle given its 
initial position and momentum, then we would fail equally miserably. What 
probably motivated the statement is the Heisenberg Inequality as applied 
to an atom of the beam. Since its transverse momentum is defined by the 
coUimation process to be a small fraction of /i/A, its position "uncertainty" 
is a large multiple of A; this is refiected in our choice of initial wave function 
ip (0; C) = 1 in the previous section, giving uniform probability for all C,. 
However, the maximum defiection, in a typical case where four even lines 
are visible either side of the central line and the transverse momentum at 
entrance to the beam is zero, is 2.4.10~^rad, so, for a laser of width 0.1mm, 
the maximum change in the value of C, at exit is 24nm or 0.04A. Although the 
observed variable is the momentum, which means that in QM parlance, is 
a "hidden" variable, I assert that it is not unreasonable to maintain that, to 
within the atomic diameter, somewhat less than Inm, each atom in the beam 
has a well defined C, which varies only slightly as the atom crosses the laser. 
For the moment I shall confine the model to the case 7 — >^ 0, which means 
we are assuming the upper internal state of the atom is infinitely short lived 
and only even lines of the momentum spectrum are seen. 

The model I propose is that of a Markov process on the set of integers 
n (r; the transition matrix being A^n (O^ that is the probability of a tran- 
sition from n to m in an interval 5t is \mn (C) <5t + o (5r) . The Markov 
property means we assume that transitions in successive intervals occur in- 
dependently. I shall make two further assumptions: [i) the process is single- 
step, so Xmn = unless m = n ± 1] this is suggested by the property (fTBj) of 
the QM process (li) the process is homogeneous, and therefore [7j additive, 
so Xm+r,n+r = A^n- With thcsc assumptious, the transition matrix has only 
two independent components, denoted Xn+i,n = « (C) and Xn-i,n = P {()■ 

To summarize, the stochastic model of the process associates, with a 
specific location ( in the standing wave, a specific Markov process with the 
parameters a {() , P {(). 
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The characteristic function for n (r; () is 

/s {9; r; C) = {e"'') = exp [ra {() {e^' - l) + r/3 (C) (e"^^ - l)] , (22) 
and its occupation probabihties are 

n/2 



exp[-aT-/?T]4(2TV^) . 



(23) 



The predicted hne intensities are obtained by integrating over the "hidden" 
variable that is 

pI ir) = lf'Pn (r; C) dC , (24) 

where P„ (r; C) = -Pn ("T"; C + ■ I shall assume further that a{( + 7r/2) = 
P (C), which results in 

P„(T;C) = P-n(r;C + ^/2) , (25) 

and hence 

Plir)^p'-nir)^l r'Pn{r;OdC. (26) 
A concrete example of such a model is 

l + sin2C . s 1 - sin 2C . . 

« (0 = ^ , P (C) = ^ ■ (27) 

This model is depicted in Fig. 4, where the directions and magnitudes of the 
transition rates a{Q ,P {() are indicated for a few different locations of the 
atom within the standing wave. Putting ^ — — tt/A, the intensities become 

2e-- 



Plir) ^ tan'^e7„(Tsin20rf^ 

Jo 



e-r ^ r (r + 1/2) r (n + r + 1/2) 
TT ^ r! (n + r)! (n + 2r)! ' ^ ^ 

and the characteristic function is (compare eqn(7)) 

oo 

Fs {e- t)^ Pn (^) = exp [-T (1 - COS ^)] Jo (r sin e) . (29) 
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Figure 4: The Kapitza-Dirac effect according to the stochastic modeL For a 
given location of the atom in the stationary wave, there is a pair of transition 
probabihties, for transverse impulses of ±2/i/A respectively. The sum of these 
probabilities is the same for all locations. For example, the two probabilities 
are equal at a node or an antinode, while one of them is zero at a point 
midway between a node and an antinode; the atom moves preferentially 
towards the nearest node. 

4 Comparison of the models 

The latter model has a variance 

{n\-\r^ + r, (30) 

as compared with the QM variance of Whilst the variances become 

indistinguishable for large r, for small r there is an essential difference; the 
initial variance is of order in the QM model, and of order r in the stochastic 
model. I postpone discussion of this disagreement to the Discussion section. 

In making a more detailed examination of the spectrum we begin by 
comparing the asymptotics of the QM and stochastic models in the limit r 
+00. The QM intensities were obtained in Section 2, and wc now compare 
them with the asymptotics of the stochastic model, which are obtained from 
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its characteristic function 



TT 



7r/2 



7r/2 



exp 



:i + sin2C) + -^-Ml-sin2C) 



namely 



Pl{r) 



1 

2wi 



Gs (r;z)z-^~^dz , 



dC, (31) 



(32) 



c 



where C is a closed contour enclosing the origin. This set of functions has the 
surprisingly simple asymptotic behaviour (obtained by selecting a steepest- 
descent contour for C) 

1 



r) 



{\n\ < t) 



(33) 



which, on averaging over the velocities of the atomic beam, gives exactly the 
same asymptotics as the QM model, that is (see Fig. 2) 



rn 



r) 



pliT')HiT,T')dT'r.pTiT) 



(34) 



Now, turning to small values of r, 1 plot, in Fig. 5, the unsmoothed 
intensities, that is and p^, of the first six lines as functions of r. Note 
that the QM and the stochastic models agree as to their orders of magnitude; 
in particular the latter model gives just four visible lines on either side of the 
centre line for the case r = 3. However, these curves show the disagreement 
of variances noted above; it shows up as a zero slope at the origin for all p^°, 
compared with a negative slope for and a positive slope for pf; the zero 
slope for higher is a consequence of the single-step assumption which we 
made in constructing the stochastic model. 

A more serious disagreement is the oscillatory dependence of pn on r in 
the QM model. Indeed that model predicts zero intensity for p„ (r) whenever 
Jn (t) has a zero. An averaging over r, as above, will give a smoothed 
intensity which never completely vanishes, but nevertheless, for small r, the 
oscillations persist even with such smoothing. On the other hand, in our 
concrete stochastic model po decreases monotonically, while the other p„ rise 
to a single maximum and then decrease monotonically, that is there are no 
zeros. This behaviour is common to the whole family of Markov models, as 
may be seen by considering the derivative 

a r 



dr 



-aT—j3T 



[a + (3) In- 2y/a(3I, 



-n+l 



(35) 
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Figure 5: The momentum spectrum of the Kapitza-Dirac effect according 
to the stochastic model. The continuous hues represent the intensity of the 
nth hne as a function of the time r spent in the laser, and the dashed line 
represents the equivalent intensity in the QM model. At r = 3 only the lines 
n < 4 are visible. 

The right hand side is negative for all positive a, /?, r and all nonnegative n, 
and therefore, substituting in (j^ . 



Thus pg (r) certainly decreases monotonically, as does also pg (r) but not 

Po (''") Po {'^) ■ For n > Q the implication of the inequality is somewhat 
more complicated, but it is certainly not satisfied by p^ (r). 

5 Inclusion of odd-momentum states 

I shall now improve the model by including the odd states, so that n takes 
half-integral as well as integral values. A jump of ±1/2 from an even (that 
is integral n) state occurs with probability 5r, that is either direction is 
equally probable, and a jump of +1/2 from an odd state occurs with prob- 
ability 5t (1 + sin2(^) /7^, while a jump of —1/2 from an odd state occurs 




(n > 0,r > 0) . 



(36a) 
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with probability 5t (1 — sin2(^) /7^. Adopting the standard classification of 
stochastic processes, the new model may be described as a pair of coupled 
additive processes, an additive process being one which is homogeneous 
and Markov. Because of the factors of 7"^, this model gives a correction to 
Fso even to zero order, but we shall see that such corrections are substantial 
only for r of order 7^. Outside of this range, the corrections to Fso are of 
order 7^ only, so they do not significantly reduce the disagreements we have 
just found between and p^. 

The characteristic function for the new model is 



- h {0; r; C) dC , (37) 

J-n/2 



where the function fs [6; r; Q may be decomposed into parts coming from 
even and odd states, that is 

00 00 
fs = fi + f2,fi = Y.Pn (r; C) e"'-' , /2 = E (r; () e^("+^/^)^ , (38) 

— CXD —00 

the P„and Pn+1/2 being the occupation probabilities of the even and odd 
states. Then /i, /2 have the time derivatives 

A = 27-2 [cos (6/2) + t sin (6/2) sin 2C] /2 - 2f, , 

A = 2 cos (^^/2) A -27-2/2 (39) 

with the initial values /i {9; 0; () = 1, /2 {0; 0; () = 0. The solution is 

(71 - 2) - (72 - 2) e-^i^ 
7i - 72 

f2 = 2cosr ^ (40) 

\2J 71-72 

where 

1 + 72 ± ^1 + 272 (cos ^ + i sin ^ sin 2C) + 7^ , 

7i,2 = 5 ■ (41) 

T 

To order 7^, and for 7^ ^ r ^ we may discard the terms in e~'^'^^ to 
obtain 



/s = /i + /2 = e" 



■720''" 



1 + ^ ( 2 cos ^ - 2 + (2r + 1) 720 - tj^o 



(42) 
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where 

720 = 1 ~ COS 9 — isinO sin 2( . 
Averaged over ( this gives 



(43) 



Fs{e-T)=Fso{9;T) 



1 + 7^ ( cos - - 1 



7' (2r + 1) 



^80 [i^]r) 



2 



The spectrum, for 7^ ^ r ^ 7 ^, is then 



(44) 



p„ = p° (1-7') - 



72 (2r + 1) .0 72r ., 



2 -Pn - —Pn . Pn+1/2 = y (Pn + Pn+l) , (45) 



and it may be extended to the range r <^ 7 ^ by adding the effect of the 
terms in e~'^^'^, that is 

2 22 
Apo = y , Api/2 = -y , ^pi = ^ e'^^" , Ap„ = (n > 1) , 

(46) 

the quantity tq having been defined in (|7I). In Fig. 6 I have plotted Pq (r) in 
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Figure 6: The intensity of the centre hne as a function of the transit time 
in the interval < r < 1. The upper curve depicts the QM model, and the 
lower curve the stochastic model. The parameter 7 has been taken as 0.2. 



the range < r < 1, including a plot of (r) for comparison. 
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The new stochastic model indicates the role of the states of odd momen- 
tum. They have an intensity of order 7^, because the transition time from 
an upper to a lower atomic state is smaller than that of the reverse transi- 
tion by a factor of order 7^. This indicates a crucial role for the zeropoint 
electromagnetic field (ZPEF) which has also played an important role in 
the program, developed by Emilio Santos and myself [H] [IHl , and designed to 
achieve a local realist understanding of the optical Bell experiments. When 
the detuning frequency exceeds the resonant Rabi frequency, "spontaneous" 
transitions, that is transitions induced by the ZPEF, are more frequent than 
laser-induced ones. Note that the new stochastic model differs radically from 
the QM model for the case that the incoming atom is in its upper state, be- 
cause in that case the fast transition, with probability proportional to 7"^, 
occurs first. This gives a spectrum with strong lines at n = ±1/2, ±3/2 . . . , 
and weak lines at n = 0, ±1, ±2 . . . , that is a reversal of the pattern shown 
for an incoming atom in its lower state. The QM model predicts no difference 
between these two spectra, so the discrepancy may provide an experimental 
method for determining which is the more correct out of the two models. 

6 Discussion 

Before discussing the disagreements between the QM and stochastic models, 
I emphasize the agreement we obtained in the asymptotic limit t —>■ 00. 
It is easily shown that the choice of a {() and P {() made in (jTfj) is, up 
to a phase shift in (, the only one which gives asymptotic agreement with 
the QM model. There is a simple explanation for this, namely that, in the 
stochastic model, the drift in an interval r is r [a {() — (3 (C)], which, with the 
choice we have made, becomes rsin2(^. This, without diffusion, is precisely 
the deterministic model occurring in Section 2 as the classical limit of the 
QM model (see Fig. 2). Hence the deterministic parts of both the QM and 
stochastic models give the same results. 

The disagreement between the models, which we found for very small r, 
may well be irrelevant, since the QM model has very rapid oscillations (see 
Fig. 3), and we have just seen that the modified Markov model also produces 
a radical change in the intensities for very small r. The quantum mechani- 
cal behaviour, whereby the initial probability of transition from a pure state 
changes from 1 by a quantity proportional to r^, is a general characteristic, 
called the Quantum Zeno Paradox, according to which a continuously ob- 
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served system cannot undergo a transition. This paradox has never received 
a satisfactory resolution. On the other hand (n^) is proportional initially to 
r for any Markov process. Note, however, that, although Fig.6, like the first 
diagram in Fig. 5, does indeed show an initial decrease in proportional to 

r, compared with for both and p^, the initial curvature of the lat- 
ter is very large compared with that of the former, which means that direct 
observation of the Zeno phenomenon would be extremely difficult. 

I turn finally to the disagreement shown in Fig.5, in particular the os- 
cillatory behaviour of P^°(t). We need to know the extent to which this 
behaviour of the line intensity is actually supported by experiment, in partic- 
ular whether the observed spectrum is consistent with ()36ap . If the existence 
of zero-intensity lines for certain r is confirmed by experimental evidence, 
then the class of stochastic processes may have to be extended to allow for 
the possibility that the atom has a memory of a recent transition having 
occurred. 

The comparison I have made between the QM and stochastic models, or 
between the wavelike "atom" of Fig.l and the more recognizably atomic ob- 
ject of Fig. 4, indicates to me that the atom of Democritus, or of Boltzmann, 
is by no means dead. Inequality (|36ap provides us with a means to determine 
which of these simple models gives the better agreement with experiment. It 
would also be interesting to try repeating the experiment with an incoming 
beam of atoms in the upper state, to see whether the pattern of strong and 
weak lines is actually inverted, as predicted in the stochastic model. 
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